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Age constraints and fine tuning in variable-mass particle models 
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VAMP (variable-mass particle) scenarios, in which the mass of the cold dark matter particles 
is a function of the scalar field responsible for the present acceleration of the Universe, have been 
proposed as a solution to the cosmic coincidence problem, since in the attractor regime both dark 
energy and dark matter scale in the same way. We find that only a narrow region in parameter 
space leads to models with viable values for the Hubble constant and dark energy density today. 
In the allowed region, the dark energy density starts to dominate around the present epoch and 
consequently such models cannot solve the coincidence problem. We show that the age of the 
Universe in this scenario is considerably higher than the age for noncoupled dark energy models, 
and conclude that more precise independent measurements of the age of the Universe would be 
useful in distinguishing between coupled and noncoupled dark energy models. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite ,1,] has confirmed that the Universe is very nearly 
flat, and that there is some form of dark energy (DE) 
that is the current dominant energy component, account- 
ing for approximately 70% of the critical density. DE 
is smoothly distributed throughout the Universe and is 
causing its present acceleration l2j.y]. It is generally mod- 
eled using a scalar field, the so-called quintessence mod- 
els, either slowly rolling towards the minimum of the po- 
tential or already trapped in this minimum 0, 0, 0-0 
WMAP has also given new observational evidence that 
nonbaryonic cold dark matter (DM), a fluid with vanish- 
ing pressure, contributes around 25% to the total energy 
density of the Universe.in agreement with a set of cos- 
mological observations 0, ul • 

An intriguing possibility is that DM particles could in- 
teract with the DE field, resulting in a time-dependent 
mass. In this scenario, dubbed VAMPs (variable-mass 
particles) the mass of the DM particles evolves ac- 
cording to some function of the dark energy field d>, s uch 
as, for example, a linear function of the field [81 H(i lllj 
with an inverse power law dark e nerg y potential or an 
exponential function [Tl hH El 111 ITrf with an exponen- 
tial dark energy potential. In this work we are particu- 
larly interested in the exponential case, since it presents 
a tracker solution; 1 that is, there is a stable attractor 
regime where the effective equation of state of DE mim- 
ics the effective equation of state of DM j^, . 

This type of solution also appears when the DE field 
with an exponential potential is not coupled to the other 
fluids '9,ll3- In fact, Liddle and Scherrer showed that for 



a noncoupled DE, the exponential potential is the only 
one that presents stable tracker solutions 0|. In this 
case, however, it is not able to explain the current acceler- 
ation of the Universe and other observational constraints, 
unless we assume that the field has not yet reached the 
fixed point regime [Tsj . 

The tracker behavior is interesting because once the 
attractor is reached, the ratio between DM energy den- 
sity p x and DE energy density p$ remains constant af- 
terwards 2 . This behavior could solve the "cosmic coinci- 
dence problem," that is: why are the DE and DM energy 
densities similar today? 

In this article, we want to study the extent of fine 
tuning required in this class of models given the ob- 
servational constraints on the Hubble constant Hq and 
the dark energy density today f^o- We also compute 
the resulting age of the Universe for the allowed mod- 
els. In Sec. ^ we present the main equations of the 
VAMP scenario and the constraints that have already 
been discussed elsewhere. In Sec. IHII we show that based 
on very robust constraints such models are not able to 
solve the cosmic coincidence problem. Section [W] dis- 
cusses the constraints imposed by the age of the Universe 
in these models, showing that better model-independent 
measurements of a lower limit would provide a possible 
way to distinguish between the noncoupled and coupled 
dark energy models. In Sec. E]we present our conclu- 
sions. 

II. THE IDEA OF EXPONENTIAL VAMPS 

In the exponential VAMP model, the potential of the 
DE scalar field <f> is given by 



V(<j>) = V e w/m *, 



(1) 



'Electronic address: urbanoOift.unesp. br 
^Electronic address: rosenfel@ift.unesp.br 

1 We adopt the nomenclature "scaling" solution to denote solu- 
tions where p^ oc a~ m and p m <x a~ n , whereas a "tracker" 
solution has m = n. 



2 This behavior does not occur in the situation in which the mass 
varies linearly with the field, although one can let the ratio be- 
tween p x and p^ remain constant for a lo ng time, until the DE 
finally comes to dominate completely llll . 
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where Vq and (3 are positive constants and m p = 
M p /V&k = 2.436 x 10 18 GeV is the reduced Planck mass 
in natural units, % = c = 1. Dark matter is modeled by 
a scalar particle \ of mass 



M x = M- 



-\(4>—(f>o)/m p 



(2) 



where M x o is the current mass of the dark matter particle 
(hereafter the index denotes the present epoch, except 
for the potential constant Vq) and A is a positive constant. 
Since dark matter must be stable, its number density n x 
obeys the equation 



h x + 3Hn x = , 



(3) 



where the dot denotes a derivative with respect to time, 
H = a/a, and a(t) is the scale factor. Assuming that the 
dark matter is nonrelativistic, its energy density is given 
by Px — ^x n xi ano - ^ follows that 



where we have assumed that the Universe is flat, 



9 = — - — 9j 

PcO 



9 



xo 



9 



bo 



9 r Q — 1 



(12) 



where p c0 = 3m 2 Hl = 8.1 h 2 x 10~ 47 GeV 4 is the present 
critical density and the Hubble parameter is Hq = lOO/i 
km s _1 Mpc -1 . pi, and p r denote the energy densities 
of baryons and radiation, respectively, which for i = r,b 
satisfy the equations pi + 3Hpi(l +uJi) — 0, where u>b = 
and u> r = 1/3. It is useful to express Vq in units of 



Pr 



4.2 x 10~ 47 GeV 4 



(13) 



the critical density value for h = 0.72. Observational 
limits on the Hubble parameter gives h — 0.72 ± 0.08 
[20|. which we adopt in what follows. 

The Friedmann and DE motion equations become, in 
terms of u, 



Px + 3Hp 



(4) 



Since the total energy-momentum tensor has to be con- 
served, the fluid equation for dark energy is given by 



H 



2 ill 



1 



3m 2 



:(Pb + Px) + Pr + 3V 



X Px (3V 



in,, 



(14) 



p + 3Hp^(l + uj^) = — p x , 



(5) 



where — p$/ ' p^ — (^<jr — V)/(\4> + V) is the usual 
equation of state parameter for an homogeneous scalar 
field. 

These equations can also be written in the form 



where 



p x + 3Hp x (l+oj^) =0 , 
p^ + SHp^l+J^) =0 , 



x 3Hm v 3m K ' 



(6) 
(7) 

(8) 



H z = 



(l/3mp (p x + Pb + p r + V) 



1- (l/6m|) <f>' 2 



(15) 



Using the fact that the right-hand side of Eq. I|14|) is 
the derivative with respect to the field 4> °f an effective 
potential [Tl| . 



U e// (0)=U(0)+p x (0) , 



(16) 



one can show that there is a fixed point value for the 
field, given by dV e f / (</>) / defy = 0: 



In 



(X + /3) (A + j8) \\p x0 e x ^/ m p 



PV 



(17) 



w * ~ ^77 =u$- (9) 

3Hm p p^ 3m p p^ 



are the effective equation of state parameters for dark 
matter and dark energy, respectively. Primes denote 
derivatives with respect to u — ln(a) = — ln(l + z), where 
z is the redshift, and ao — 1. From the above equations 
one can also obtain the DE equation of motion, 



3H<j>- 



X Px 



pv 

m p 



(10) 



The Friedmann equation for a Universe with DE, DM, 
baryons, and radiation is given by 



1 



3m p 



Pb + Pr + P X + \4> 2 + V((f>) 



(11) 



At the present epoch the energy density of the Universe is 
divided essentially between dark energy and dark matter. 
In this limit, using the above solution, one obtains 



9(1-. — 1 — 9^ 



A 



(A + /3) 2 A + /3 



(e) 



A + /3 



(18) 



(19) 



which is a stable attractor for /3 > -A/2 + (\/A 2 + 12)/2 
|12lll6| . The equality between lu x and uj^ in the attractor 
regime comes from the tracker behavior of the exponen- 
tial potential 0,0,0] in this regime. 

The density parameters for the components of the Uni- 
verse and the effective equations of state for the DE and 
DM for a typical solution are shown in Fig. ^ Notice 
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FIG. 1: Top panel: Density parameters of the components of 
the Universe as a function of u — — ln(l + z) for A = 3, (3 = 2, 
and Vo = 0.1p c . After a transient period of baryonic matter 
domination (dot-dashed line) , DE comes to dominate and the 
ratio between the DE (solid line) and DM (dashed line) energy 
densities remains constant. Bottom panel: Effective equations 
of state for DE (solid line) and DM (dashed line) for the same 
parameters used in top panel. In the tracker regime both 
equations of state are negative. 



that the transition to the tracker behavior in this exam- 
ple is currently occurring. 

Let us assume observational upper and lower limits to 
the DE equation of state parameters uj l obs < lu 1 ^ < uj^ bs 
from SNIa, for instance. Then, from Eq. (|19|l it follows 
a constraint on the parameter space (A, /3) given by 



/?<(>) 



1 



Mi), 

J obs I 



U(I) | 

J obs I 



(20) 



which is valid if the field has already reached the attractor 
solution. 

Analogously, in the attractor solution, the density pa- 
rameter fi^o can also constrain the parameter space, since 
Eq. ((TH1) implies 



(3 > (<) li ' A 



2 n 



(DE)u(l) 
obs 



(21) 



a»[(2 n£*w> 



-i) -wT Ml) 



-1)] + 12 



(DE)u(l) 



(DE)u(l) 
obs 



where f2 
value of DE density parameter. 



2 n 



is the upper (lower) limit to the current 
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FIG. 2: Constraints on parameter space (A, (3) given by DE 
equation of state and DE density parameter. The upper 
(lower) limit on equation of state allows the regions under 
(above) the dot-dashed (dotted) line for u) b s = —0.4 (—0.8). 
On the other hand, the upper (lower) limit on the DE density 



parameter, Vt 



(DE) 
obs 



0.8 (0.6), requires the parameters to be 



above (under) the solid (dashed) line. The hatched region is 
allowed by both constraints. 



Amendola [T^ has studied in detail the observational 
constraints imposed in this scenario by the dark energy 
equation of state (obtained from luminosity distance con- 
sistent with high z supernovae) and density parameters 
when the attractor regime has already been reached. He 
obtained —0.8 < < —0.4 at 95% confidence limit 
(C.L.). Figure |3 shows the allowed region in the (A,/3) 
parameter space using this result together with a conser- 
vative bound 0.6 < tt<f,o < 0.8 on the DE density param- 
eter. 



III. FINE TUNING 

Tracker solutions of noncoupled exponential 
quintessence models in which the dark energy equation 
of state tracks that of the background are not able to 
accelerate the Universe today. Exponential VAMPs 
solve this problem because the background acquires a 
negative equation of state. However, it is necessary to 
verify whether fine tuning of the model is necessary to 
explain the observations today. 

In order to quantify the amount of fine tuning required, 
we have solved numerically the coupled equations of mo- 
tion for the scalar field, varying Vo, A, and /3 in the region 
A = [0.01,20], = [0.01,20], and V = [Q.01/o c , 1.5p c ], 
with step sizes AA = A/3 = 0.25 and AV = 0.05p c , 
generating about 1.9 x 10 5 models, from which we have 
selected those that agree with observed cosmological pa- 
rameters. 

Figure |3 shows the normalized number of models that 
satisfy some observational constraints as function of Vo, 
given in units of p c . The dashed line shows the nor- 
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FIG. 3: Normalized number of models that satisfy the ob- 
servational constraints discussed in the text. 
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FIG. 4: Density parameters of the components of the Uni- 
verse as function of u = ln(l + z) for A = 3, j3 = 2. In the top 
panel Vo = W~ 4 p c , and in the bottom panel Vo = 10 3 p^. 



malized number of models that satisf y t he Hubble pa- 
rameter constraint, h = 0.72 ± 0.08 [2(J]. Notice that 
in order for this constraint to be satisfied it requires 
Vq to be of the order of p c , the critical density today. 
The dot-dashed line is obtained requiring both the Hub- 
ble parameter constraint and the equation of state con- 
straint, -0.8 < w$ < -0.4 at 95% C.L. Q. The 
dotted line corresponds to the case when the Hubble 
parameter constraint and parameter density constraint, 
0.6 < fl^Q < 0.8, are taked into account, and the solid 
line when these three parameters are constrained simul- 
taneously. One can see that the narrow allowed range 
of V = [0.25p c ,0.45p c ], 68% C.L. (F = [0.15p c ,0.55p c ], 
95% C.L.), is essentially determined by the constraints 
given by h and Sl^o, which are very conservative and ro- 
bust. Hence, fine tuning in this class of models seems 
unavoidable. 

What does this fine tuning in Vq mean? Vo is the main 
parameter that determines the epoch in which the field 
reaches the attractor regime. As can be seen in Fig. ^ 
for the allowed values of Vo, the Universe is entering the 
attractor regime around the present epoch. This implies 
that the equation of state and the parameter density are 
still varying today, although their current values are very 
similar to the values at the fixed point. 

However, if one varies Vq a similar behavior of CIde 
and CIdm sets in earlier or later, as can be seen in Fig. 
QJ However, unrealistic values of the Hubble parameter 
result in this case. For instance, for Vb = lO 3 ^, where 
the fixed point is reached earlier, we obtain h ss 10. In 
this sense, the cosmic coincidence problem is not solved 
in these models, since one has to choose the overall scale 
of the potential in order to obtain realistic cosmological 
solutions, very much like noncoupled dark energy models. 



IV. AGE OF THE UNIVERSE 

The potential to use the age of the Universe as a con- 
straint to dark energy models has been discussed in sev- 
eral recent papers [2ll HI HI El • The WMAP collabo- 
ration pj has given a stringent value for the age of the 
Universe, to = 13.7 ±0.2 Gyr at 68% C.L. However, this 
result is model dependent and obtained from direct inte- 
gration of the Friedmann equation for the running spec- 
tral index ACDM model (cold dark matter model with 
a cosmological constant). Observations also indicate a 
lower limit for the age of the oldest globular clusters (and 
consequently for the age of the Universe) of to = 10.4 Gyr 
at 95% C.L. H2. 

Using the limits discussed above and the position of 
the first Doppler peak in the WMAP experiment it is 
possible to put an upper limit on the equation of state 
(assumed constant) of the dark energy, < —0.67 at 
90% C.L. |23. However, current upper limits on the age 
of the Universe alone cannot constrain the lower value 
of the equation of state, since for large values of the age 
(to > 18 Gyr) the cosmology is essentially independent of 
the value of the equation of state [23| , unless we assume 
that reionization took place too early in the Universe 
history (as indicated by WMAP data), and consequently 
the upper limit for the age of the Universe is close to the 
upper limits coming from globular clusters, to ~ 16 Gyr, 
as in the case discussed in Ref. |22^ . 

A model independent approach to estimating the age 
of the Universe has been performed by means of a 
parametrization for the equation of state of dark energy 
that allows it to vary [24| . Using cosmic microwave back- 
ground (CMB) and supernova data the age of the Uni- 
verse was found to be to — 13.8 ± 0.3 Gyr at 68% C. L. 
Even this model independent approach, however, is not 
able to parametrize the VAMP scenario, since in this case 
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FIG. 5: Constraints os parameter space (A, 0) imposed by the 
age of the Universe for Vo = 0.1p c , Vo = 0-3p c , and Vo = 0-5p c . 
Lines correspond to the same limits from Fig. [5] 
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FIG. 6: Age of the Universe versus the present dark energy 
equation of state for models that satisfy the Hubble parameter 
constraint. 



the effective equation of state of the CDM is also variable, 
and the strong limits in the age of the Universe obtained 
by these CMB analyses cannot be used to constrain the 
models discussed here. For this reason, we have used only 
the conservative limits obtained from globular clusters to 
study the potential of the age of the Universe to constrain 
this scenario. Therefore we adopt 10.7 < to < 16.3 Gyr 
at 95% C.L. H3. We have added to the age of globu- 
lar clusters 0.3 Gyr, since the WMAP data indicates an 
early reionization of the Universe P, . The age of the 
Universe in the VAMP scenario is given by 



to = H 1 J \ f^ exp 



+ n M (i + zf + fi ro(l + *) 

r 



[i+j:\u)]du 



(22) 



O x o exp 



i+4 e) (") 



du 



-1/2 



dz 



(l + z) 



where itj is the initial value for u (ui = —30 for our 
numerical calculations) and Zi its corresponding redshift. 

To explore the region of parameter space that is able to 
satisfy the age constraint we have used the same proce- 
dure described in the preceding section, varying A and 
in the region [0.01, 20], with a step size AA = A/3 = 0.05 
and a fixed value of Vo, generating about 1.6 x 10 5 mod- 
els. Models were required to satisfy the constraints on 
the Hubble parameter and on the age of the Universe. 
Results are shown in Fig. for three different values of 
Vo. 

The constraints imposed by the age of the Universe 
are in the same direction of those imposed by the dark 
energy density and the equation of state, indicating a cor- 
relation among them. This can be seen explicitly in Fig. 
HO Notice that such a correlation is exactly the opposite 
that occurs when the equation of sta te ( for noncoupled 
dark energy) is taken to be constant [22|: higher values 



(e) 

of (jJ^q imply lower ages of the Universe. 

In the constant case, the relation between the age and 
ui is direct: for more negative equations of state the con- 
tribution of dark energy becomes important later (since 
p scales as a~ 3 ( 1+ ^). Consequently, for a longer time 
the Universe is matter dominated, and its age is lower 
compared with a Universe in which the dark energy has 
a more positive equation of state. 

For the VAMP scenario, however, the situation is a lit- 
tle more complicated: models that satisfy the age of the 
Universe consist of a positive correlation in the param- 
eter space that has a lower slope for larger Vo (Fig. EJl. 
For larger Vo, the Universe reaches the attractor regime 
earlier, and consequently has a larger value of the age 
of the Universe (for the same j3 and Hq), since the Uni- 
verse is accelerating for a longer time. Besides that, for 
the same [3, larger values of Vo imply larger values of A 
for the models that satisfy the age constraint, and conse- 
quently, from Eq. (|19fl , the equation of state tends to — 1 . 
Despite this correlation, due the large degeneracy in the 
allowed region, it seems to be very difficult to constrain 

uj^q based on limits on to and vice versa. 

Figure [7| shows the age of the Universe for the models 
that satisfy the Hubble parameter and the dark energy 
density observational constraints. Here we have used step 
sizes AA = A/3 = 0.2 and AVq = 0.05p c for the region 
A = [0.01, 20], /3 = [0.01, 20], V = [0.1p e , 0.8p c ]. 

Fitting the distribution of models as function of ages, 
the age of the Universe in this VAMP scenario was found 
to be 



to = 15.3+^ Gyr at 68% C.L. 



(23) 



which is considerably higher than the age of models of 
noncoupled dark energy |l],[24|. This seems natural, since 
in these models the CDM also has an effective negative 
equation of state and accelerates the Universe. Thus, 
measurements of the age of the Universe could help to dis- 
tinguish between coupled and noncoupled models. This 
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FIG. 7: Distribution of models that satisfy Hubble and DE 
density parameters constraints as a function of age of the 
Universe. Fit to the points gives to = 15.3^'^- Gyr at 68% 
C.L. 

result is very conservative, since it relies only on the well 
established limits of the Hubble constant and the dark 
energy density today. 

One can see that lower limits on the age of the Uni- 
verse are not useful to constraint these models. Upper 
limits on the age, however, are potentially interesting. 
As an example, we can speculate what would change if 
the limits on the age of globular clusters pll were 
symmetric: to = 12-5 ± 2.2 Gyr, and add to this age 0.3 
Gyr for the formation of these objects. In this case, the 
upper limit on the age of the Universe would be 15 Gyr, 
which would exclude a large part of the models allowed 
by the present limits. The parameter space for these lim- 
its (Fig. |HJ is much more constrained than the previous 
case, practically excluding Vq > 0.5p c . 

V. CONCLUSION 

The VAMP scenario is attractive since it could solve 
the problems of exponential dark energy, giving rise to a 
solution of the cosmic coincidence problem. However, in 
order to obtain solutions that can provide realistic cosmo- 
logical parameters, the constant Vb has to be extremely 
fine tuned in the range Vb = [0.25/5^, 0.45p c ] at 68% C.L.. 
This implies that the attractor is being reached around 



the present epoch. In this sense, the model is not able to 
solve the coincidence problem. 

We have found that there is a negative correlation be- 
tween the value of the age of the Universe today and the 
current equation of state of dark energy: larger values 

of oj^q correspond to a lower age of the Universe. This 
result is the opposite of the case in which the equation 
of state of dark energy is constant. 

A generic feature of this class of models is that the Uni- 
verse is older than noncoupled dark energy models. In 
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FIG. 8: Constraints os parameter space (A, /3) imposed by 
to = 12.5 ± 2.2 Gyr, and adding 0.3 Gyr as the age of the 
Universe in the epoch the globular cluster were formed, for 
Vo = O.lp^, Vb = 0.3p c , and Vq — 0.5p c . Lines correspond to 
the same limits from Fig. [5] 



fact, we found to = 15.3 jlo'7 Gyr. Better model indepen- 
dent determination of the age of the Universe could help 
to distinguish among different contenders for explaining 
the origin of the dark energy. 
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